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Abstract 

We present the results of a detailed study of the large transverse momentum 
Drell-Yan process, pp — > (7*, Z*)X l^l~X at collider energies, with either one 
or both protons polarised, allowing the study of single- and double-spin asymmetries 
respectively. We show how these asymmetries obtained from angular distributions 
of the leptons in the 7* (or Z*) rest frame, can be used to get information on 
the polarised parton distributions. Numerical results for the asymmetries and the 
cross-sections are presented, and the sensitivity of the asymmetries to the initial 
parton distributions indicates that these can be used as effective probes of the spin 
structure of the proton. 
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1. Introduction 



Over the last two decades, deep-inelastic scattering experiments have given us important 
clues to the structure of the nucleon. Many of the predictions of the QCD-improved 
parton model have now been tested in these experiments, and we have a reasonable 
understanding of the structure of the nucleon in terms of its constituents - quarks and 
gluons. This understanding, however, is far from complete. There are several aspects 
of hadronic phenomenology that do not yield to a complete description in terms of the 
QCD-improved parton model. One striking example of this is the description of the spin 
structure of the nucleon. 

The limited amount of experimental information on the spin structure of the nucleon 
that is available is from polarised deep-inelastic scattering experiments - scattering of 
polarised leptons off polarised nucleon targets. Over five years ago, the EMC collaboration 
first published its results ^ obtained from polarised muon-proton scattering experiments. 
For the polarised structure function gi{x,Q'^), which the EMC experiment measured in 
the range 0.015 < a; < 0.7, the first moment gives 

Tl = f\f{x,Q^)dx = 0.12S± 0.032. (1) 
Jo 

This integral was obtained by assuming a smooth extrapolation based on Regge behaviour 
for the unmeasured low-x region x < 0.015. In the parton model, we have the relation 

9f{x) = W^"Mi^)- (2) 

In the above equation, Aq^{x) = (^^(a;) — qL{x)), where g£ refer to densities of quarks 
with ± helicity in a proton with helicity |, i.e. the Ag-^'s are twice the contribution to 
the nucleon's spin of a quark of flavour /. Data from hyperon decays determine two 
independent linear combinations of the first moments of Ag-^'s. Using the information 
in Eq. |l| with that available from hyperon decays, it is possible to determine the first 
moments of Ag-^'s {i = u, d, s) separately. This yields for the integral of twice the sum of 
the spins of the quarks 

AS = y A^dx = J{Au + Ad + As)dx = 0.12 ± 0.17. (3) 

This is the famous EMC result that the total spin carried by the quarks is small and is 
actually compatible with zero. 

More recently, data from the SMC muon-deuteron scattering experiment 0] and from 
the SLAC muon-^He scattering experiment |^ have become available. Taken at face value, 
the SLAC results seem to be in disagreement with the EMC and SMC results. However, 
a careful analysis of all the uncertainties show that the values of AS obtained from 
the three sets of data are not incompatible with each other. A global average yields 

AS = 0.22 ±0.10. (4) 
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This value is not very different from the value in Eq. ^; and the global average also gives 
a rather small value for AS. In the naive parton model, the theoretical problem is to 
understand the smallness of this quantity which could arise from a large, negative strange 
quark polarisation 0. However, in QCD, there is an alternative interpretation |^ of the 
data on Ti. Via the axial anomaly, there arises a gluonic contribution to gi{x) and to 

so that the conclusion that AE is small can be avoided, if the gluonic contribution is 
large enough, i.e. if the polarised gluon distribution function is non-negligible. 

It thus becomes important to have other independent information on the polarised 
quark and gluon densities. In this paper, we study the large transverse momentum Drell- 
Yan process, taking into account Z-interference, as a probe of the gluon polarisation in 
the nucleon. We use a range of polarised parton distributions, fitted to the data from 
the EMC experiment. We study this process for the energies {^/s ~ 500 GeV) planned 
at the proposed polarised p-p experiment at the Relativistic Heavy Ion Collider (RHIC) 
in Brookhaven. A detailed account of the proposal can be found in Ref. |p. At RHIC, 
experiments with both incident protons polarised and those with only one of the incident 
protons polarised, are envisaged. The process that we are interested in is pp + X, 

where the lepton pair has large transverse momentum and where either one or both initial 
protons are polarised, allowing the study of single-spin and double-spin asymmetries, 
respectively. For this process, we study the angular distributions of the leptons, and 
demonstrate how this is sensitive to the polarised parton distributions. At the high 
energies under consideration, virtual photons and virtual Z's contribute to lepton pair 
production. The parity-non-conserving Z vertex is crucial to our considerations because 
this results in a non-zero single-spin asymmetry. 

Some of the results discussed in this paper, namely those relating to the simplest 



asymmetries of the cross- sect ions, were published in an earlier letter |]T0|. In this paper 



we discuss in detail the theoretical framework for obtaining the angular distributions, and 
present complete expressions for the additional asymmetries associated with the multipole 



parameters. Numerical results beyond that presented in Ref. [10|, are also included in the 
present paper. The remainder of this paper is organised as follows : in Sec. 2, we begin 
with a discussion of the explanation of the EMC experiment in terms of the axial anomaly 
and a large gluon polarisation and emphasise the need to have independent experimental 
information on the magnitude and the shape of the polarised gluon distributions. Af- 
ter briefly discussing the proposals to measure the gluon polarisation that exist in the 
literature, we discuss the conditions under which a non-vanishing single-spin asymmetry 
can be obtained. In Sec. 3, we discuss the process pp — > + X and give the general 
framework for obtaining the angular distributions of the leptons in the rest frame of the 
vector boson (7* or Z*), and in Sec. 4 we present the numerical results for the asymme- 
tries and the cross-sections. We state our conclusions in Sec. 5. In Appendix A, we list 
the complete expressions for the density matrices and the multipole parameters, and give 
the expressions for the decay parameters in Appendix B. 
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2. Gluon polarisation and the single-spin asymmetry 



We begin this section with a brief description of the anomaly explanation ||^ of polarised 
deep-inelastic scattering experiments. In the naive parton model, the flavour- singlet part 
of the first moment of the polarised structure functions is given by 

^llsinglet = ^ (5) 

However, in the QCD-improved parton model, this relation is no longer true: there is 
an additional piece due to the interaction of the photon with a gluon via a quark loop. 
Naively, this contribution would be expected to be small, as it involves an extra factor of 
as as compared with the lowest-order photon-quark diagram. But the contribution is not 
small because the longitudinally polarised gluon distribution evolves in such a way as to 
compensate exactly the logarithmic decrease with of as- For the first moment of the 
singlet part of gf, the gluon contribution is given by the triangle axial anomaly, so the 
singlet part of the measured polarised structure function includes a gluon contribution, 
which is given by the following relation 

nUngiet = ^AS - ^NjAg, (6) 

where Nf is the number of quark flavours that circulate in the loop. The empirical 
smallness of Ti\singiet is explained in this picture by postulating a cancellation between 
the AE and the Ag contributions; in order that this cancellation occurs one requires a 
large gluon polarisation. 

It is important to realise that the anomaly explanation of the EMC experiment does 
not provide any estimate of Ag, neither do we have any other theoretical reason to believe 
that it must be large. An independent experimental determination of this quantity is 
crucial, and will, in fact, help us understand the spin structure of the nucleon. There 



have been several suggestions in the literature |TI| on how the gluon polarisation may be 
experimentally determined. Most of these are aimed at experiments where the beam and 
the target are both longitudinally polarised - the asymmetries that can be studied in this 
case are the double-spin asymmetries. From the experimental point of view, however, it 
would be much simpler if only one of the initial particles were polarised. There exist in 
the literature few suggestions on how the single-spin asymmetries thus measured may be 
used to determine the polarised parton densities. 

The most direct way of measuring a single-spin asymmetry is to measure the polari- 
sation of a final-state particle, and the measured asymmetry would then simply correlate 
the final polarisation to the initial polarisation. The measurement of the polarisation of 
a final-state photon in direct photon production or of a final-state lepton in Drell- 



Yan dilepton production [T^ have, therefore, been suggested. However, it is not an easy 



task to measure the polarisation of a final-state particle at high energies, and, hence, 
it becomes necessary to look for other ways of obtaining a finite single-spin asymmetry. 
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One way is to measure at least two final-state momenta and couple them through their 
cross-product to the spin of the initial particle. The cross-product and the spin are both 
axial vectors and yield a scalar correlation; this correlation, however, is T-odd and can, 
therefore, be induced only through loop corrections. A detailed study of one-loop cor- 
rections to polarised Drell-Yan has been made |T^, in order to estimate the size of the 
asymmetry arising from the imaginary part of the amplitudes. The only other way of 
obtaining a non-vanishing single-spin asymmetry is to consider parity-violating processes 
like W production |16[. This asymmetry obtains because the initial spin couples with 



the axial part of the ly-fermion vertex. In this paper, we consider another parity- violating 
process: large pt Drell-Yan with 7*-Z* interference, and study both single- and double- 
spin asymmetries. This we do by studying the angular distributions of the lepton pair in 
the vector boson rest-frame. 



3. Large-pr Drell-Yan with ^*—Z* interference 

The subprocesses that contribute to the Drell-Yan process at large pr are qq annihi- 
lation and qg Compton scattering: 

q + q Z* + g, 

q{q)+9 ^7*, Z* + q{q), (7) 

followed by 7*, Z* ^ . The Compton subprocess dominates over the annihilation 
subprocess in the case of p-p scattering. Consequently, this process is sensitive to gluon 
distributions in the p-p case. We consider polarised p-p collisions (with the first proton 
labelled A and the second labelled B). The corresponding partons from protons A and 
B are labelled a and h and these carry momentum fractions Xi and X2, respectively. We 
will consider the case where both A and B are polarised and also the case where A is 
polarised and B is not. 

It is useful to think of the Drell-Yan process in terms of a production process where 
the photon or the Z is produced, and a "decay" process for the boson into an pair. 
We analyse the decay angular distributions of the lepton in the rest frame of the 7* (or 
the Z*) with the z-axis taken to be the direction of the momentum of the 7* or the Z* 
(in the cm. frame). In this frame, 9i and 0; are the polar and azimuthal angles of the 
lepton. 

We can factorise the Feynman amplitude for the Drell-Yan process pp l^l^X into 
two parts: one part of the amplitude specifying the production of the virtual vector boson 
(7* or Z*), and the other describing the decay of this boson into a pair of leptons. Thus, 
we may write (with /~, standing for the helicities of the lepton and the antilepton, 
respectively) 



M^' = — 



MA^;p.(7)Ai+'l7) + x{M)M>^,,p^{Z)A'^'- {Z) 



(8) 
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where the Ma^;po-(7,Z) are the production amphtudes and A'-^'- {^,Z) the decay amph- 
tudes for the 7* and the Z*, and A, /x, p, a are the hehcities of the vector- meson, the final 
parton and the initial partons h and a, respectively; M is the invariant mass of the lepton 
pair and the function x{^) is given as 



X{M) 



M2 



M2-M|+ iMTz' 



(9) 



where Mz and Tz are the mass and the width of the Z, respectively. Squaring the 
amplitude in Eq. (|]) and summing over all helicities except that of the vector boson, we 
arrive at the following expression : 



E, 



1 



+rpllA{^'~{l)A^y'^*{Z) + XPlAMr'-(Z)Ai:'-*(7) 



7^t Al + l~ 



where the unnormalised density matrices pw are defined as 
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E/.''(^i)/i'(^2)M,,^,.(7)M* (7) 



(10) 



:iii 



etc. The factors /^(xi), fp{x2) are the appropriate density functions for partons of 
helicity p. We write the Feynman amplitude for the decay as 



A 



i+i- 



no sum on a) 



:i2) 



where « = /" — /"*" = ±1 only, because for fast leptons only the helicities (H ) or ( — h) 

are allowed. The M^s are essentially standard model coupling constants. Then expanding 
the density matrices in terms of multipole parameters, and integrating over 0; yields 



CXh) + \x\'C^ti*{Z) + 2Re[x*Cftl*i^Z)] 



1=0 



Yioiei), (13) 



where Yiq are the spherical harmonics and the Ci are decay parameters [|13, and we have 



t: 



I* 



E(lA|lA;/0)p 



AA, 



Ci = V3j2{lO\la-l-a)Ml. 



(14) 



In the above equation stands for |Mq,(7)P and \Ma{Z)\'^ for the pure 7 and Z terms 
and Ma{'y)M*{Z) for the interference terms. 
Finally Eq. (|13|) can be rewritten in the form 



\M\ 



i) 



M2 



(15) 



1=0 
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where the Di^s are: 



Do = tM + {v' + a')\x\XiZ) + 2v{Rex*K{^Z), 
D, = -V6a[v\x\HliZ) + iRex*)tlilZ)], 

D2 = -^[tl{^) + {v' + a')\x\Hl{Z) + 2v{Rex*)tl{7Z)]. (16) 

Here v and a are the vector and the axial couphngs of the Z to the leptons. The expressions 
for the multipole parameters for the annihilation and the Compton subprocesses, which 
involve the parton distributions and the couplings, are presented in Appendix A. The 
decay parameters C; are discussed in Appendix B. 

4. Numerical results 

For the Drell-Yan process AB — > I'^l'CX, (C is an associated quark- or gluon-initiated 
jet) where the pair has transverse momentum pt and rapidity j/i, the differential cross- 
section to produce the lepton at angle 9i (in the reference frame described in the previous 
section) is then given as 



da 



AB-^l+l-CX 



-= / dx^^TY^DiYi^iOi), (17) 
I ■' 1=0 



dpTdyidrdcosO, 

where r = M^/s, with s the cm. energy, and Xi the momentum fraction of proton A 
carried by parton a. The factor JF is given by 

T= ^ (18) 



xi - ^^x^ + Arey^ 



where xt = IptI \f^- 

For each setting of the spins of the coUiding protons the differential cross-section in 
Qi is controlled by the independent parameters -Do,i,2, each of which contains information 
about the parton distributions. To isolate information on the polarised distributions, it is 
necessary to form asymmetries, either by reversing the spin direction of one of the protons 
with the other being unpolarised (single-spin asymmetries) , or reversing the spin direction 
of one of the protons, the other being polarised (double-spin asymmetries). 

For each spin setting, the £)j's can be projected out from a knowledge of the angular 
distribution : 

^ nz J AB^l+l-CX 

Di^^ , , , , . YM)smeidei, (19) 

Clearly Dq is the simplest, being essentially the cross-section to produce an pair with 
given Pt, yi and r. 
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The full set of asymmetries that one can, in principle, measure is the following : 

Dq — Dq da^ — d<7_ 

As 



+ Z^o da+ + da- 



where the labels +, — refer to the helicities of the protons A and B, and for j = 1,2 

Ai 



Dt+ - D-+ 



- 4+ +4+- ^^^^ 

In what follows, we shall present results for all these asymmetries. We shall assume that 
the jet recoiling from the lepton pair is not detected and so we must integrate over Xi. 
We illustrate our results for the case yi = and fixed values of M, and show how the 
asymmetries vary with px- We take ^/s = 500 GeV to correspond to RHIC energies. 

In order to study the sensitivity of our asymmetries to the polarised parton distribu- 
tions we utilise the following range of models : 



Set I : Ag large. As = 0, 

Set II : Ag and As both moderately large. 

Set III : A^i = 0, As large. 



which were described in Ref. [|T8| and which all fit the EMC data on gi{x). For the 
unpolarised distributions we use Owens' Set 1.1 distributions [|19|. 

In Fig. 1, we show the calculated asymmetries for M = 10 and 50 GeV and for 
M = Mz and with 61 integrated over. The figures in the top row are predictions for 
single-spin asymmetries, the bottom row for double-spin asymmetries. In each figure, the 
three curves correspond to the three different sets of polarised densities used. We see that 
for small dilepton masses the single-spin asymmetries are very small, but that they are 
appreciable for large dilepton masses. This is because the single-spin asymmetry arises 
from the parity-violating nature of the Z-fermion coupling, and hence is large only when 
the effects of Z- interference are important. The double-spin asymmetries are large over 
the whole range of dilepton masses. In Fig. 1, there is a strong dependence upon the 
choice of polarised parton distributions and the asymmetries, if measured with reasonable 
accuracy, will be very helpful in teaching us to what extent the various partons in a proton 
are polarised. In Fig. 2 we have plotted the cross-section as a function of pt for the three 
values of dilepton mass given above. In Figs. 3 and 4 we have plotted the asymmetries 

d U = l!2). Since the structure of the contributions to D2 is very similar to those 
of -Do contribution, the corresponding asymmetries also look very similar in shape. The 
single-spin asymmetries A\''^ are large only close to M = Mz- 
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5. Conclusions 



Data on the polarised structure function gi{x) in deep-inelastic scattering, from the EMC, 
SMC and SLAC experiments, if interpreted in the naive parton model, suggest that only 
a small fraction of the proton's spin arises from the spin of its quark constituents. One 
explanation has the contribution of the valence quarks cancelled by an unexpectedly large 
strange quark contribution. Another approach argues that the naive parton model for- 
mulae are invalidated by the axial anomaly and explains the data in terms of a large 
gluon contribution to the proton's spin. Independent information on the magnitude of 
the strange quark and gluon polarised distributions is thus of great importance and needed 
urgently. In this paper we have studied both the single- and the double- spin asymmetries 
of the parameters which describe the angular distribution of the lepton in the large trans- 
verse momentum Drell-Yan process at RHIC energies (v^ = 500 GeV). At these high 
energies, Z*'s can contribute significantly to the Drell-Yan process, inducing a parity- 
violating single-spin asymmetry. This asymmetry may be measured in experiments where 
only one of the initial particles is polarised. 

We have discussed the framework in which we compute the angular distributions of 
the leptons in the vector-boson rest-frame, and have presented the complete results for 
the multipole parameters, which describe the decay angular distributions of the lepton. 
Using these angular distributions, we have constructed single- and double-spin asymme- 
tries, which we then studied numerically using three sets of polarised parton distributions 
that are consistent with the deep-inelastic scattering data. We find that the single-spin 
asymmetries are measurably large only at large dilepton masses, i.e. close to the Z peak. 
The double-spin asymmetries are large even at smaller dilepton masses. Both the single- 
spin and the double-spin asymmetries are sensitive to the polarised gluon distribution and 
may be used as probes of the spin structure of the proton. 
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Appendix A 



In this appendix, we present the expressions for the density matrices and the multipole 
parameters, for the annihilation and the Compton processes. 

The couphngs of the 7* to both quarks and leptons is given by the interaction term 
—ieQfj^, while the coupling of the Z* to the fermions is given by the term ie'-y^ivj — aj'-)^), 
where as usual 



-'3 



(22) 

In the above equations. Is is the third component of the weak isospin, 6w the Weinberg 
angle and Q/, Qj are the charges in units of e. In what follows, we use v and a to denote 
the Z* couphngs to leptons and Vf and a/ to quarks of flavour /. The results are presented 
in terms of the variables f and 6, where 9 is the subprocess scattering angle, and 

T 



X1X2 



(23) 



The unnormalised density matrices for Z* production via the annihilation process 
qq Z*g are given as : 



pn(^) = Ka{l + f2) Y.{M^M^2) - Aqf{x^)Aqf{x2)) [{vj + a}) (cot 



tan 



,9' 



+2vfaf I cot - — tan - 



,9 



+ {Mf{xi)qf{x2) - qf{xi)Aqf{x2)) 



,9' 



.9 



9' 



Poo{Z) 



Pi-i{Z) 



{y^ + a^) I cot - — tan - 1 + 2vfaf ( cot^- + tan 

^a8fE|(g/(a;i)g/(a;2) - Aqf{xi)Aqf{x2)){v) + a]) 
f ^ 



+{Aqf(xi)qf(x2) - qf(xi)Aqf(x2))(2vfaf)j, 

^a4f E|(?/(^i)?/(^2) - Aqf{xi)Aqf{x2)){v} + a}) 
f ^ 



Pw{Z) 



+ {Aqf{xi)qf{x2) - qf{xi)Aqf{x2)){2vfaf)j, 

KaV2^{l + f) ^|(g/(a;i)g/(a;2) - Aqf{xi)Aqf{x2)) {vj + a}) ^tan^ - cot 

- {^Qf{xi)qf{x2) - qf{xi)Aqf{x2)) x 



(9 9 
—2vfaf I tan — h cot- 



9 



[Vf + ttf) I tan- + cot- 1 - 2vfaf I tan- - cot- 1 |, 



Po-iiZ) = KaV2f(l + f)J2{iQf(xi)(lf{^2)-Mf{xi)Aqf{x2))\{vj + a}) ( cot^ - tan^ 



9 9' 



—2vfaf ^cot- + tan^ 



- {^(If{xi)qf{x2) - qf{xi)Aqf{x2)) 



Vj: + a^) I tan- + cot- 1 + /w/a/ I tan- — cot 



= Ka{l + T^) ^ {qf{x^)qf{x2) - Aqf{xi)Aqf{x2)) 
f ^ 

+2vfaf [ tan - — cot - 



[vj + a}) (tan^^ + cot^^ 



+ {^<lf{xi)qfix2) - qf{xi)Aqf{x2)) 

2 2\ / 2^ 2^ \ f-> / 2^ 2^ \ 1 

•jj^ + (^f) I tan - — cot 2 j ^"^Z*^/ I tan - + cot - I |, 

where the sum runs over all flavours. The overall factor Ka is given by 

4 eV 



(24) 



is:, 



" 9(l-f)2- 



(25) 



The unnormalised density matrices for Z* production via the Compton process qg ^ Z*q 
are given as: 



Pu{Z) 



K,^\G{x,)qf{x2) 
f 



(vf - ajf ni - f)2sec^^ + f^cos^^tan^^ j + [vj + affcos^^- 



+G{xi)Aqf{x2) 
+AG{xi)qf{x2) 
+AG{xi)Aqf{x2) 



[Vf - a/) I (1 - r) sec - + r cos -tan ^j-{vf + af) cos - 

/ \9//-i '~\9 9^ ^^9 9^ 4.^1 / \9 9^ 

{vf - a/) I (1 - ^) sec - - T cos -tan ^j+{vf + af) cos - 



(vf - aff (1 - ff 



sec 



-^2 2^, 4^ 

r cos -tan 



poo(^) = K,2fsm^^-J2{iG{xi)qf{x2)-AG{xi)Aqf{x2))2{vj + a}) 



-Avfaf{G{xi)Aqf{x2) - AG{xi)qf{x2))j, 
p,_,{Z) = K,Hin^^Y.{i-G{xi)qf{x2)+AG{xi)Aqf{x2))2{vj + a}) 
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+Avfaf{G{xi)Aqf{x2) - AG{xi)qf{x2)) \, 



e . e 



(vf + a/)2] + {AG{xi)qf{x2) - G{x,)Aqf{x2)) (^{vf - aff fUn^^- - {vf + aff^ |, 
XeV2f cos^sin^ Y,{{G{x^)qf{x2) - AG{x,)Aqf{x2)) ({vf - aff + {vf + ajf 

ftan^^) + {G{x{)Aqf{x2) - AG{x^)qf{x2)) {{vf - aff - {vf + affrtan^^^ |, 



P-i-i{Z) 



K,^\G{x,)qf{x2) 
f 



+G{xi)Aqf{x2) 
+AG(xi)qf(x2) 
+AG{xi)Aqf{x2) 



e. J 



e 



-{vf + o/) (1 - t) W- + f^cos^-tan^- ] + {vf - afYcos'- 



.e' 
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-{'^f + «/) (1 - ^) sec - - f cos -tan -]- {vf - af) cos - 



{vf + QfY ( (1 - rfsec'^- - f"cos"-tan^- 1 - (^;/ - a/)^cos"- 



(26) 



where, again the sum is over all flavours, and the overall factor Kc is given as 

1 e^g^ 



6(l-f)' 



(27) 



We can construct the unnormalised multipole parameters using the expressions for the 
density matrices given above. For the annihilation process, we obtain 

f,{Z) = 2K,J2{M^M^2) - ^1fi^i)^1fi^2)){vj + a}) + {Aqf{x,)qf{x2) - qf{x,)Aqf{x2)) 
f ^ 



2vfaf 



(l + f^) ftan^^ + cot2^ j +4f 



tl{Z) = V2Ka J2{ {'lf{xi)Qf{^2) - Aqf{xi)Aqf{x2))2vfaf + {Aqf{xi)qf{x2) - qf{xi)Aqf{x2)) 
f 



{vj + a}) 



(1 + f^) Ut^^-tan^^ 



tliZ) 



2\/]^^aE{(^/(^i)9/(^2) - Aqf{x,)Aqf{x2mvj + a}) + 
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{Aqf{xi)qf{x2) - qf{xi)Aqf{x2))2vfaf 
and for the Compton process, we get 



(1 + f2) (tan^^ + cot^^ 



(28) 



tliz) 



ir,^2(t;J + 4)jG(xi)g/(x2) 
/ 



(1 - ffsec^- + cos^- 1 + fhan'^- + 2fsm^- 



+AG{xi)Aqf{x2) 



(1 - f)^sec^- - cos - 1 + retail - 



2fsin - 



^vjafi G{xi)Aqf{x2) 
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(1 - f)^sec - + cos - 1 + f^tan"^- + 2fsm - 



+AG{xi)qf{x2) 
1 

+AG{xi)Aqf{x2) 



e 



(1 — f) sec^ cos^- I 1 + f^tan 

^ ^ 2 2 \ 2 



2fsin^ 



,9 



-j=Kc -Avfaf^G{xi)qf{x2) 



9 



,9 



(1 - f)^sec - + cos - retail''- - 1 



x 2 2 ^ 2 ^ 

— r sec cos - 

^ 2 2 



+2{vj + aj)\Gixi)Aqf{x2) 



fhan^- - 1 



9 



.9 



(1 - f)^sec^- + cos^- f^tan*- - 1 



(1 _ ff^eA - cos^- { f^tan^- - l' 
^ ' 2 2 I 2 , 



+AG{xi)qf{x2) 

^ K,Y.2{v) + a])\G{x,)qj{x2) 
f 



(1 - 'r)^sec^- + cos^- f 1 + f^tan'^- 
^ ' 2 2 V 2, 



4fsin^ 



9' 



+AG{xi)Aqf{x2) 



,9 



9 



(1 - f)^sec^ cos - 1 + retail''- + 4fsin^- 



4:Vfafl G{xi)Aqf{x2) 



,9 



9 



(1 - f)^sec - + cos - 1 + f^tan"^- - 4fsin - 



+AG{xi)qf(x2) 



,9 



9 



,9 



(1 - f)^sec'' cos - 1 + retail''- + 4fsin - 



(29) 



The multipole parameters given are only for quark-gluon Compton scattering. The cor- 
responding multipole parameters for antiquark-gluon Compton scattering are given by 



-af 



(30) 



in Eq. 

So far, we have hsted only the Z* production density matrices and multipole param- 
eters. We need to have similar expressions for 7* productions and the expressions for the 
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interference terms. These expressions can be easily obtained from the above expressions 
for Z* by simple replacements of the couplings, which we summarise in the following 
rules : 

1. To obtain the 7* multipole parameters, the replacements a/ — > 0, f/ Q/, should 
be made in the expressions for the Z* production multipole parameters. 

2. To obtain the multipole parameters corresponding to the interference terms, the 
replacements vj+aj — > —2QfVf, ^Vfttf —Qfttf, should be made in the expressions 
for the Z* production multipole parameters. 

Having listed all the multipole parameters, we are now in a position to write down the 
complete expressions for the D/s defined in Eq. To present these in a compact form, 
we define the following combinations of the couplings 

af = Q}-2vvfQfRex* + {v^ + a^){vj + a})\x\\ 

Pf = 2afVf[{v' + a')\x\'-vQfRex*], 

7/ = aa/[2w/|xP - Q/Rex*], 

= v{vj + aj)\x\'^ -VfQfRex*, 



(31) 



l-fl^sec^-i 



2 2 ^ 2 ^ 
1 — f) sec - ± cos - 
^ 2 2 

.9 



and the functions, 

R±{f,o) = 
s±{r,e) = 

Uif,9) = 

Vif,9) = 
We may now write each Dj as 



9^99 
cos — h rsm-tan- 
2 2 2 



1 2 



r tan 



,9 



:i+f^ 



:i + f^ 



2 ^ 2 ^ 

tan — h cot - 
2 2 

2 ^ 2 ^ 
tan cot - 



+ 4f. 



= E[^f + + Df] + [9 
f 



TT 



01 



(32) 



(33) 



where the sum is over flavours, and by — tt — ^] is intended to imply the following 
symmetrization 

fA{xi)gB{x2)F{9) ^ gA{x^)fB{x2)F{'K - 9), (34) 
where /, g are parton densities pre-multiplying some function F of 9. 
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Then we have, for the annihilation subprocess, 

Df = 2KaU{f,e)iy{qf{xi)qf{x2) - Aqf{xi)Aqf{x2))af + {Aqf{xi)qf{x2) - qf{xi)Aqf{x2 
Df = 2V3KaV{T, 0)l^{qf{xi)qf{x2) - Aqf{x,)Aqf{x2)hf + {Aqf{x^)qf{x2) - qf{xi)Aqf{ 

Df = y|x„[[/(f,^)-12f]{(g7(xi)g/(a:2)-Ag/(xi)Ag^(x2))a/ 

+ (Ag/(xi)g/(x2) - g/(xi)Agy(x2))/3/|, 
and for the Compton subprocess, 

Df = 2K,i^[G(x,)qf(x2)R+(rJ) + ^G(xi)Aqf(x2)R-(rJ)]»f 
-[G{xi)Aqf{x2)R+{f,9) + AG{xi)qf{x2)R-{f,9)]f3fy 

Df = 2V^K,^^[G{x,)qf{x2)S+{T,e) + AG{x^)Aqf{x2)S4T,e)]^f 
-[G{xi)Aqf{x2)S+{f, e) + AG'(xi)g/(a;2)-S_(f, 

Df = ^i^e{[G(:ri)g/(:^^2)r+(f,^) + AG(a:i)Agy(a;2)r_(f,^)]a/ 

-[G(xi)Ag;(x2)T+(f, 0^) + AG{x,)qf{x2)T_{T, (36) 

Df^s are obtained from Df^s by making the replacements /9/ — > — /3/ and 7/ — > —7/. 

Appendix B 

In this appendix, we discuss the decay process V where V denotes 7* or Z*. 

We first list the Feynman amphtudes A^^^ , where l^, l~ refer to the helicities of the 
antilepton and the lepton, respectively, and A to the helicity of the vector-boson. Since 
the vector-boson is massive (of mass M), we compute these amplitudes in the vector- 
boson rest frame. In this frame, we take the polar and azimuthal angles of the lepton to 
be 61 and 0/, respectively. 

For fast leptons, because of the vector and the axial- vector couplings, only the helicities 
(-1-,-) or (-,-|-) can occur and the Feynman amplitude has the form 

4^'- = M^e^f^'dlM (37) 

where a — l~ — l'^ — ±1. 
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Upon calculating the Feynman amplitudes, we find that 

M„ = V2iMp{a) (38) 

where l3{a) depends only on the electroweak coupling constants. For the decay of the 
photon one finds 

(3^^\l)=ie (39) 

and, by parity invariance, we get 

/3W(-1) =/3W(l) (40) 
Then, from Eq. |14|, the decay parameters Ci are given by 

C;^)(00; 00) = V3{2M^) J2 \f3^''\a)WlO\la; 1 - a); (41) 

a 

which gives 

Cj^)(00;00) = 4e2M2, 
C{^^(00;00) = 0, 

C^^)(00;00) = 2726^ Ml (42) 
For the Z*, /3^^\l) = ie{v — a) and f3'^^\—l) = ie{v + a), and so we get 

Cf)(00;00) = {v^ + a^)C^^\00;0Q), 



Cf^(00;00) = -AVqM^-^ 



e va. 



Cf)(00;00) = {v^ + a^)Ci^\00;00). (43) 

Finally, for the interference term, we obtain 

Ci^^^(00;00) = vC^^\00;00), 
Ci^^)(00;00) = -2VQM^e\ 

C5^^^(00;00) = vCi''\00;00). (44) 
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Figure captions 



Fig. 1 Asymmetries for the Drell-Yan process at y/s = 500 GeV, as a function of pt- The 
asymmetries shown in this figure are constructed out of the cross sections differential 
in Pt, yi, and f (the jet rapidities have been integrated over), with f fixed at M/ y/s 
and Hi = 0. The cross sections are integrated over the lepton angle 6i. The dilepton 
mass M is shown on top of the figures. The full, dashed and dotted lines correspond 
respectively to Set I, II and III polarised densities, described in the text. The upper 
row of plots are for the single-spin asymmetries As, the lower for the double-spin 
asymmetries Ad- Note that the scale for the first two plots is much smaller than for 
the rest. 

Fig. 2 Cross-sections for the Drell-Yan process at = 500 GeV, as a function of pr- The 
caption da/dT refers to da / dpTdyidf , with f fixed at M/y^ and yi = 0. The cross 
sections are integrated over the lepton angle 9i. The full, dashed and dotted lines 
correspond respectively to M=10 GeV, M=50 GeV and M = Mz, where M is the 
dilepton mass. The cross sections are in nb/GeV. 

Fig. 3 Same as in Fig. 1, but for the asymmetries Al and A]^ defined in Eq. 

Fig. 4 Same as in Fig. 1, but for the asymmetries Aj. and defined in Eq. 
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